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The scale of Baryon Acoustic Oscillations (BAO) imprinted in the matter power spectrum provides 
an almost-perfect standard ruler: it only suffers sub-percent deviations from fixed comoving length 
due to non-linear effects. We study the BAO shift in the large Horndeski class of gravitational 
theories and compute its magnitude in momentum space using second-order perturbation theory 
and a peak-background split. The standard prediction is affected by the modified linear growth, as 
well as by non-linear gravitational effects that alter the mode-coupling kernel. For covariant Galileon 
models, we find a 14-45% enhancement of the BAO shift with respect to standard gravity and a 
distinct time evolution depending on the parameters. Despite the larger values, the shift remains well 
below the forecasted precision of next-generation galaxy surveys. Models that produce significant 
BAO shift would cause large redshift-space distortions or affect the bispectrum considerably. Our 
computation therefore validates the use of the BAO scale as a comoving standard ruler for tests of 
general dark energy models. 
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I. MOTIVATION 

One of the most exciting promises of modern 
cosmology is the possibility of testing fundamental 
physics using the largest scales available to observa¬ 
tion |T]. Among other developments, the signatures 
of Baryon Acoustic Oscillations (BAO) have pro¬ 
vided an invaluable test of models for cosmic acceler¬ 
ation through their imprint in the Cosmic Microwave 
Background [2] and the distribution of Large Scale 
Structure (LSS) j3 using either galaxies @H8], the 
Lyman-cr forest or quasars El 113 (see Ref. |TH - 
m for reviews). To an excellent approximation, the 
BAO signal in the LSS provides a comoving stan¬ 
dard ruler that traces the expansion of the universe 
and probes the onset of cosmic acceleration. 

Non-linear corrections are known to introduce a 
small departure from the perfect standard ruler be¬ 
havior, systematically shifting the BAO scale to¬ 
wards smaller values at low redshift. This effect has 
been well studied for cold dark matter cosmologies 
with a cosmological constant using perturbation the¬ 
ory P3H22 and simulations [1221 f2T| (for earlier works 
see (25, 2Bj). The result is that the BAO scale im¬ 
printed in the matter distribution shrinks by approx¬ 
imately 0.3% at redshift zero [22} [24]. However, this 
value relies on the assumption that gravity is New¬ 
tonian in the scales of interest. 


Little attention has been devoted to the non-linear 
BAO evolution in more general theories of gravity. 
Since the shift in the BAO scale is comparable to 
the sub-percent level of precision expected by forth¬ 
coming galaxy surveys m that aim to test such the¬ 
ories, it will be necessary to understand the effects 
of non-standard gravity on the BAO scale to cor¬ 
rectly interpret the data in the next generation of 
dark energy experiments. 

II. PEAK-BACKGROUND SPLIT 
COMPUTATION OF THE BAO SHIFT 

Sherwin and Zaldarriaga have explained the BAO 
shift in terms of the effect of long modes on the 
short scale power spectrum [ 21] . In their picture, 
large overdense regions undergo less overall expan¬ 
sion, reducing the size of the physical BAO scale 
with respect to the average (see also [2B]). This effect 
is not compensated by underdense regions, because 
cosmic structures in overdense regions undergo more 
growth and give a larger contribution to the power 
spectrum. Therefore, local differences in expansion 
and growth lead to a net shortening of the comoving 
BAO scale, causing a small departure from the stan¬ 
dard ruler behavior. Alternativelly, the shift of the 
BAO scale can also be understood as arising from 
contributions to the power spectrum which are off- 


phase with respect to the linear prediction mm- 
The BAO shift can be studied by comparing the 
non-linear power spectrum to a rescaled version of 
the linear one [15]: 

P(k) « Pu(k/a) = Pu(fc) - (a - 1 )kP[ 1 {k) -\ -, 

(1) 

where the shift can be read from the coefficient 
of the second term (more sophisticated templates 
are actually used to obtain the BAO scale from 
data, but we will stick to this description for sim¬ 
plicity). Here and below the power spectrum is 
defined as {S(k)S(k')) = (2ir) 3 5D{k + k')P(k) and 
Pnm oc {S n 5 m ), where we expand the density con¬ 
trast as 5 = 8 i + 62 + ■ ■ ■ ■ 

One can compare Eq. with the prediction from 
standard perturbation theory 


P(k) = (P11 • • • + Pm) + (P22 • • • + Pmu) ■ ( 2 ) 

All P\ n contributions are proportional to Pn(fc) and 
thus do not contribute to the second term in Eq. 0 
m- Only the mode-coupling terms (second paren¬ 
thesis) do contribute to the shift, with the first of 
such contributions given by 


P 22 (fc) 


d 3 q 

(2P)3 


P 2 (fc 



Pu(k-q)Pii(q) ■ 


(3) 

Here P 2 is the second-order symmetrized mode¬ 
coupling kernel [29] . The rest of the computation in 
the peak-background split approximation proceeds 
by expanding in k/q, integrating with a cutoff at 
kBAO & n d extracting the coefficient of kP[ 1 (k) from 
the result (see |2T] for further details). The long 
modes with q <C k ~ ksAO describe the effect of the 
large fluctuations on the smaller scales. 

The computation can be generalized to alternative 
theories of gravity by noting that the structure of the 
kernel F 2 is preserved on sub-horizon scales, but each 
term acquires a time-dependent coefficient Ci(t) 


P2(p,q) = Co + Civ (J + p) + ° 2 - l) ’ W 

which reduces to the standard constant values, Cq = 
17/21, C\ = 1/2 and C 2 = 2/7, under matter domi¬ 
nation in the case of standard gravity (we drop the 
time-dependence for notation convenience). Explicit 
computations in the sub-horizon, quasi-static limit 
of Horndeski theories with non-relativistic matter 
determine that the modifications to the kernel co¬ 
efficients are not independent EDI El] and satisfy 

Ci = \, C 0 + 2 -C 2 = l. (5) 


The BAO shift can be read by plugging the gen¬ 
eralized kernel Q, into the mode-coupling power 
spectrum |3]), expanding to leading order in q/k, per¬ 
forming the integration and comparing with Eq. 0 
This generalizes the result of Ref. eh to 

a - 1 = ^ (2C0 - ^ (Si) , ( 6 ) 

where Eqs. ([5]) have been used to write the result in 
terms of the monopole Co- In the above expression 
the integration over the momentum leads to the long 
mode variance 

rk B AO 

( 6 l)= 7— ^47r g 2 Pn(g,f) « <j 2 rs {t ), (7) 

Jo \ Z7T ) 

where we use a cut-off at BAO scale, well estimated 
by the sound horizon at the drag epoch kgAO ~ 
1 /r s ( Zd ) E2] ■ Following Ref. [21] , we use the square 
of the variance of the density field on a sphere of 
radius r s (zd) for the computation of the BAO shift: 
<*i> « < 7 / (t). This gives a slight underestimate 
with respect to the shift measured in simulations 
of standard cosmology, but we expect comparison 
among models to be accurate. 

LSS surveys observe galaxies, which are known 
to be biased with respect to the underlying matter 
distribution. The effects of non-linear density-halo 
bias can be parameterized as 

Sh(x) = M(x) + ^b 2 S 2 (x) + • • • . (8) 

where the bias parameters bi, b 2 relate the mat¬ 
ter and the halo overdensity (the above expansion 
should hold on large scales). The halo-halo power 
spectrum generalizing Eq. ([2| reads 

P h (k) = b\ (Pn + P 22 ) + bib 2 P sj 52 + • • • , (9) 

where second-order terms that do not involve mode 
coupling have been omitted as they do not con¬ 
tribute to the BAO shift. The first parenthesis con¬ 
tains the linear and mode-coupling matter power 
spectrum renormalized by the linear bias. The sec¬ 
ond term mixes second order corrections with non¬ 
linear bias, and affects the BAO shift. This term is 
given by 

p 6 2 1 8 2 (k) = j ~^^2F 2 (k - q,q)Pn(k - q)Pn{q), 

( 10 ) 
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where the only difference with Eq. ^ is that the 
kernel appears linearly. As before, one can com¬ 
pute the leading order expansion in q/k and com¬ 
pare with a rescaled version of the linear halo power 
spectrum (cf. Eq. |l]) multipled by b\). Identifying 
the kP[ 1 (k) term yields the contribution to the shift, 
which now reads 



This result generalizes equation ([6]). Note that 
the relations ([5]) for the perturbation theory kernels 
make the bias contribution to the shift independent 
of the theory of gravity, i.e. of any departure in the 
value of Co- For this reason we will not consider 
non-linear bias in the next section. 


III. BAO SHIFT IN ALTERNATIVE 
THEORIES OF GRAVITY 

We focus our analysis on theories within the Horn- 
deski Lagrangian [33.,, which contains many exam¬ 
ples of interest for cosmology including Brans-Dicke, 
/(i?), Chameleons, Kinetic Gravity Braiding and co¬ 
variant Galileons. Horndeski’s theory also contains 
the characteristic interactions that appear in consis¬ 
tent theories of massive gravity and higher dimen¬ 
sional theories, and it is thus expected to effectively 
describe some of their distinctive features [31] Q 

Although our analysis is general, for the sake of 
simplicity we will present results for a covariant 
Galileon model [H] (see also |451 :15]). We fix the 
Galileon Lagrangian parameters and the cosmolog¬ 
ical parameters to the best-fit models obtained by 
Barreira et al. [T7] (without massive neutrinos), 
which have zero cosmological constant. We noticed 
that the quintic galileon model we present has a gra¬ 
dient instability in the tensor sector. However we 
decided to include it in our analysis since it is a 
good fit for the data and has interesting properties 
at second-order in perturbation theory. Indeed, as 
we shall see it produces large modifications of the 
dark matter kernel, which can be detected by study¬ 
ing the bispectrum with current surveys [48, 149| . 


1 We will not consider viable extensions of Horndeski’s theory 
[351 |36I . nor full theories containing interacting gravitons 
[33H39]. See Refs. I40ll43l for reviews on the cosmology of 
alternative theories of gravity 


Model 

(7 t> s 

2C 0 

Olk - 1 [%] 

A 

0.067 

1.62 

0.20 

Cubic 

0.071 (7%) 

1.61 (-0.4%) 

0.23 (14%) 

Quartic 

0.073 (9%) 

1.58 (-2%) 

0.23 (15%) 

Quintic 

0.071 (7%) 

1.92 (19%) 

0.29 (45%) 


Table I. Density contrast, mode coupling kernel 
monopole and BAO shift in the matter power spectrum 
for reference model and selected Galileon models at red- 
shift zero (cf. Fig. [I]). Values in parenthesis indicate the 
relative deviation with respect to a Cosmological Con¬ 
stant model. Tracer bias can be added using Eq. (111. 


Simpler scalar-tensor theories (such as Brans- 
Dicke, Chameleons or f(R)) lead to very constrained 
modifications and can not produce large contribu¬ 
tions to the Kernel m- For such models the only 
sizeable contributions to the BAO shift stem from 
the enhancement of linear growth, and would thus be 
in conflict with measurements of LSS clustering. If 
one further demands that these theories are screened 
in the Galaxy or the Solar system, the range of the 
scalar force is too short to even affect cosmological 
scales in the linear regime m- An exception is given 
by non-universal couplings to matter: Coupled dark 
matter models can significantly increase the shift of 
the BAO scale ED- 

The generalized Sherwin-Zaldarriaga formula ([6]) 
depends on the theory of gravity in two ways: a 
correction from linear physics, given by , and a 
modification of the mode coupling kernel & , given 
by Co- We compute the evolution of the background, 
the linear power spectrum and the density contrast 
of using a modified version of the CLASS code 
\Pl\ :33J based on the general description of Horn- 
deski perturbations presented in [54] (see also ED)- 
The computation of the non-linear corrections to 
the kernel follows the approach of ref. |3T] (see also 
[56] ) by taking the sub-horizon approximation and 
the quasi-static approximation (valid for covariant 
Galileons on the scales of interest |57|). We will 
also assume that the scale at which the model be¬ 
comes strongly coupled is smaller than the BAO 
scale. This is indeed the case for cubic and quar- 
tic galileon models, as suggested by a comparison 
between fully non-linear and linearized N-body sim¬ 
ulations for Galileons [581159] . 

The results for the BAO interesting quantities at 
redshift zero are presented in Table [I] All the models 
considered tend to increase the density contrast oy s 
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Figure 1. Time evolution of the BAO shift in the matter power spectrum (left panel) and integrated kernel (right 
panel) for standard and Galileon gravity models (cf. Table [I|. Red lines indicate expected sensitivity of Da across 
redshift bins from an optimistic BAO survey m (see discussion around Eq. ©)■ 


due to an enhanced effective force of gravity and 
the different background expansion. The non-linear 
corrections to Co are highly dependent on the model 
parameters, acquiring positive and negative sign and 
ranging from sub-percent in the cubic, percent in 
the quartic, and becoming fairly large in the quintic 
example. 

The time evolution of the BAO shift and the 
mode-coupling corrections are displayed in Figure 
[T] Departures with respect to standard gravity oc¬ 
cur only at low redshift and become largest in the 
accelerated era when the scalar field energy density 
drives the cosmic expansion. Besides this general 
trend, each model is characterized by a specific time 
dependence. Our results allow to distinguish be¬ 
tween a very soft non-linear regime in which the 
mode coupling is mostly determined by interactions 
of the matter fluid (standard gravity, cubic model) 
and large non-linear effects, as in the case of the 
quintic model, with the quartic case being an inter¬ 
mediate example. This is a consequence of the non¬ 
linear gravitational interactions introduced in Horn- 
deski’s theory. 


IV. DISCUSSION 

Our results show an enhancement of the BAO shift 
with respect to the standard prediction and provide 
the first estimate of this effect for modified gravity. 
It is possible to compare the predicted shift to the 
forecasted sensitivity of next-generation galaxy sur¬ 
veys. Let us focus on measurements on the BAO 


scale transverse to the line of sight Qbao = r s /DA, 
where Da is the comoving angular diameter distance 
(comparison with line-of-sight BAO yields similar re¬ 
sults). One can compare the two sources of uncer¬ 
tainty 


A P A __ + A'O 


( 12 ) 


where the first term is the observational error (as¬ 
suming known r s ) and the second term is the sys¬ 
tematic error induced by the shift, A r s /r s = a — 1. 
Weinberg et al. have provided an example forecast 
for an all-sky BAO survey in which the expected er¬ 
ror ranges from 2.8% at z = 0.15 to 0.1% 2 > 3.5 
(these data can be found in Table 2 of Ref. |12|){^] 
Figure [l] compares both terms in Eq. (12) and shows 
that the BAO shift is well below the precision for 
all the examples considered at any fiducial redshift. 
Note that the forecasted precision is mainly limited 
by survey volume, implying that more sophisticated 
observational setups will not be able to reduce the 
errors considerably. More realistic forecasts based 
on specific surveys lead to lower precision (see Ref. 

123 ). 

It is very unlikely that models more general than 
the ones considered here can lead to sufficiently large 


2 Their forecast also assumes density field reconstruction im¬ 
provements in the non-linear damping by a factor of 2. 
Since this procedure has not been validated for general the¬ 
ories of gravity, we take the forecasted precision as an op¬ 
timistic bound. 
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shifts to bias the BAO scale measurements while re¬ 
maining compatible with other observations. The 
theoretical prediction, Eq. (j6j) , allows one to iden¬ 
tify two contributions to the shift: the modified lin¬ 
ear growth and the non-linear gravitational effects 
that modify the mode-coupling kernel. Any theory 
of gravity with a very large shift requires a consid¬ 
erable enhancement of at least one of these contri¬ 
butions, which can be probed by observables other 
than BAO0 

A large departure of oy a would be ruled out by 
redshift space distortions or other clustering mea¬ 
surements. Similarly, large corrections to the mode 
coupling kernel would induce large distortions in 
the bispectrum (note that Eq. (|5j) implies that a 
> 23.5% increase in Co would change of sign in 
the quadrupole term in F 2 ). We emphasize that 
these non-linear gravitational effects are exclusive 
of fully fledged Horndeski theories (cf. quartic and 
quintic example Galileons considered here) and very 
suppressed in simpler scalar-tensor theories (e.g. 
Brans-Dicke, f(R )) or cubic theories (e.g. our cu¬ 
bic example, Kinetic Gravity Braiding m and lim¬ 
its of extra-dimensional theories EB). Most works 
on higher order perturbation theory for modified 
gravity have focused on the latter type of models 

mmrnm- 

Our findings validate the use of BAO measure¬ 
ments as a comoving standard ruler for current and 
next-generation LSS surveys, at all redshifts of inter¬ 
est and even for the most extreme theories of gravity. 
There are several refinements that can improve our 
calculation, such as including higher order perturba¬ 
tion theory corrections. Other developments would 
be necessary in order to better connect these results 
with observations, such as the inclusion of more so¬ 
phisticated bias models (which has been shown to 
affect the magnitude and time evolution |24| ~l and 
redshift-space distortions (which typically increase 
the magnitude of the shift relative to real space). 
Finally, extending our result would allow to confirm 
the validity of density field reconstruction [B5] [BB] of 
BAO for general theories of gravity]^] Despite possi- 


3 Another possibility is that a modification of gravity en¬ 
hances the BAO shift by producing a large value of the 
non-linear halo bias 62 / 61 - Such a possibility would how¬ 
ever rely on the details of halo formation and its study 
would require methods other than cosmological perturba¬ 
tion theory. 

4 The validity of BAO reconstruction schemes has been ar- 


ble refinements, the smallness of the effects ensures 
the validity of our conclusions regarding the shift. 

These are some initial steps in understanding the 
interplay between extended theories of gravity and 
the BAO scale imprints on the distribution of large 
scale structure. Further work should address other 
aspects of LSS and BAO in general theories of grav¬ 
ity in order to optimize the performance and model 
independence of the next-generation of dark energy 
experiments. This will ultimately shed light on the 
optimal strategy to test gravitational physics using 
future LSS surveys and learn more on the connec¬ 
tions between fundamental physics and cosmology. 
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